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Abstract 

We consider classical billiards on surfaces of constant curvature, 
where the charged billiard ball is exposed to a homogeneous, station- 
ary magnetic field perpendicular to the surface. 

We establish sufficient conditions for hyperbolicity of the billiard 
dynamics, and give lower estimation for the Lyapunov exponent. This 
extends our recent results for non-magnetic billiards on surfaces of 
constant curvature. Using these conditions, we construct large classes 
of magnetic billiard tables with positive Lyapunov exponents on the 
plane, on the sphere and on the hyperbolic plane. 
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1 Introduction and the Statement of Main 
Results 



In the present work we consider the billiards on two-dimensional surfaces M 
of constant Gaussian curvature r in the presence of a homogeneous magnetic 
field of magnitude (3, which is perpendicular to M. Inside the billiard domain 
Q the pointlike particle of unit charge and mass moves at unit velocity along 
curves of constant geodesic curvature (3 and reflects elastically at the bound- 
ary dQ. In the following, we will call these dynamical systems as magnetic 
billiards and M as magnetic surface if (3 ^ 0. 



Magnetic billiards on the plane have been considered in many works [BR 



H, [ BK |, [ |l'al| and on the hyperbolic plane in [ |la2| | . The study of such bil- 
liards is strongly motivated by mesoscopic physics, where such billiard mod- 
els are used as simplified version of the mesoscopic devices in the presence of 
magnetic fields. In the present paper we treat the magnetic billiards simul- 
taneously on all surfaces of constant curvature (sphere, plane and hyperbolic 
plane). For all values of r and (3 we establish a common criterion for hyper- 
bolicity of the billiard dynamics, whose geometric realization depends only 
on the type of linearized dynamics (geometric optics) on M. This extends 
our recent results ||GSG|| for non-magnetic surfaces of constant curvature. 

The dynamics on M depend crucially both on the curvature of the surface 
and on the strength of magnetic field. Firstly, let us consider the case (3 = 0. 
On the plane the neighboring trajectories separate only linearly with time, 
so that the motion of the point mass between collisions with the boundary is 
neutral. Exponential separation of billiard trajectories can only occur if the 
successive reflections from the boundaries introduce sufficient instability. On 
the hyperbolic plane the negative curvature induces exponential divergence 
of geodesies. Thus, the boundary of hyperbolic billiard can be neutral (i. 
e., with zero curvature), and the chaotic dynamics will be provided by the 
metric. On the sphere, in contrast, any two geodesies intersect twice at focal 
points. Thus, the boundary reflections have to compensate for the focusing 
effect of the sphere, in order to produce chaotic dynamics. We will call the 
mentioned above types of linearized dynamics arising on the plane, hyperbolic 
plane and sphere as parabolic , hyperbolic and elliptic respectively. 

In the presence of a constant magnetic field, an additional focusing effect 
appears. To simplify matters, let us discuss first the planar case when (3^0. 
Consider an infinitesimal beam of trajectories which emerges from the same 
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point at the time t = (fig. la). Then, by elementary calculations (see e.g., 
0) the curvature of the infinitesimal beam at the time t = s is given by 



XO) = -Pcot(ps) 



Consider now the reflection of infinitesimal beams at the boundary. Let m 
be the bouncing point, let C be the osculating circle at m and let m! be the 
second point in which the particle trajectory intersects C. Denote by X-i 
X+ the curvatures of the infinitesimal beam immediately before and after 
reflection. Then, the change of the curvature under the reflection |[ral|| , |K| 
at the bouncing point m (see fig. lb) is given by 



X+-X- = 2/3 cot (/3d), 



(2) 



where 2d = mm' is the signed length of particle trajectory between m and m! 
(when the curvature of the boundary at m is positive, 2d is simply the time 
which the particle spends after (before) reflection in the osculating circle 
C). It is a simple observation, that eqs. (1), (2) are actually the same 
ones, which have appeared in ||GSG|| for the non-magnetic billiards on the 
sphere of radius /3 _1 , where the parameters s and d have the same geometric 
meaning. Thus, the geometric optics (or linearized dynamics) on the plane 
in the presence of a magnetic field (3 is equivalent to the geometric optics 
on the non- magnetic sphere of radius (3 . More generally, we demonstrate 
in the body of the paper the equivalence of geometric optics on the surfaces 
of constant Gaussian curvature r in the presence of magnetic field f3 to the 
geometric optics on the non-magnetic surfaces of constant Gaussian curvature 
r e ff = r + p 2 . As a consequence, the type of linearized dynamics (elliptic, 
parabolic or hyperbolic) depends only on the sign of r e ff. We will call the 
parameter r e ff as effective curvature of the surface and will refer to the cases 
r e ff = 0, r e ff < 0, r e ff > as Parabolic (P), Hyperbolic (H) and Elliptic 
(E) respectively. 

Up to now, the study of hyperbolic billiards has been mainly restricted 
to the Euclidean plane (see [[Tarj|| for review). See, however, |[Ta2|| for some 



results on chaotic billiards on the hyperbolic plane, and ||Vetl||, |Vet2| , [|KSS 



for some results on hyperbolic billiards on a general Riemannian surface. In 
our recent work ||GSG|| we have generalized Wojtkowski's criterion of hyper- 
bolicity for planar billiards [ Wo2 | to billiards on arbitrary surfaces of constant 
curvature. On the basis of the equivalence of geometric optics on magnetic 



3 



and non-magnetic surfaces of constant curvature, we extend in the present 



paper the criterion of |GSGfl to the case of magnetic billiards. 

The hyperbolicity criterion in [ |GSG|| can be formulated in terms of a 
special class of trajectories, which generalize two-periodic orbits. Let Q be a 
billiard table on a surface of constant curvature. The billiard map <fi : V —>■ V 
acts on the phase space V, which consists of pairs v = (m, 9). Here m is the 
position of the ball on the boundary dQ of Q, and 9 is the angle between the 
outgoing velocity and the tangent to dQ at m. The billiard map preserves 
a natural probability measure /i on V. We denote the images of v after n 
iterations by (m n+ i,9 n+ i) = (f> n (v). The trajectory 0"(t>) is a generalized 
two-periodic trajectory (g.t.p.t.) if the following conditions are satisfied: 

1. The incidence angle and the curvature of the boundary k„ at the bouncing 
points have period 2: 9 2n = 9 2 , 9 2n +i = #i, K 2n = k 2 , k 2 „+i = «i; 

2. The length of trajectory between consecutive bouncing points is constant: 
s = m^m n+1 (see fig. 2). 

If 9i = it/2, the g.t.p.t. is an usual two-periodic orbit. 

Along a g.t.p.t. the linearized map D v (p is two-periodic, and the stability 
of a g.t.p.t. is determined by D v <p 2 . For each surface of constant curvature, 
the stability type of a g.t.p.t. is completely determined by the triplet of 
parameters (di, d 2 , s), where 2c?! (resp. 2d 2 ) is the signed length of the chord 
generated by the intersection of the trajectory m\m 2 with the osculating 
circle at m\ (resp. m 2 ). We shall use the symbol T(di,d 2 ,s) for the g.t.p.t. 
with parameters (d\, d 2 , s). 

Let us consider g.t.p.t.s for planar, non-magnetic billiards in some detail. 
Here s is the euclidean distance between consecutive bouncing points, and 
di = Ti sin 9i, i = 1,2, where are the radii of curvature of the boundary dQ 
at the respective points. If the curvature of the boundary at the bouncing 
point is zero we take r« = — oo as the radius of curvature and di = — oo 
respectively. By an elementary computation, T{d\,d 2 ,s) is unstable if and 
only if 

{[di, d 2 ] U [d\ + d 2 , oo) if d\, d 2 > 
[0,oo) if di, d 2 < (3) 

[0,di + d 2 ] U [di,oo) ifdi>0,d 2 <0. 

Moreover, the trajectory is hyperbolic (i. e., strictly unstable) if s is in the 
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interior of the corresponding interval, and the trajectory is parabolic if s 
is a boundary point (in the limiting case di = d 2 = — oo the trajectory is 
parabolic for any value of s). 

We consider two classes of unstable g.t.p.t.s. The g.t.p.t. T(di,d 2 ,s) is 
B-unstable if in eq. (3) s belongs to a "big interval" : 

![di + d 2 , oo) if di, di > 
[0,oo) ifdi,rf 2 <0 (4) 

[di,oo) if di > 0, d 2 < 0. 

On the contrary, if s belongs to a "small interval", then T(di,d 2 ,s) is S- 
unstable: 

J [di,d 2 ] if di,d 2 > , , 

1 [0, di + d 2 ] if d 1 >0,d 2 < 0. 1 J 

Note that a small interval shrinks to a point when | d\ \ = \d 2 \. 



It has been demonstrated in [PSGf , that the notions of B-unstable and S- 



unstable g.t.p.t.s are generalized to arbitrary surfaces of constant curvature, 
where the analogs of (3,4,5) exist. The concept of g.t.p.t.s and the associated 
structures make sense for billiard on any surface immersed in a magnetic field. 
Since the stability properties of the trajectories depend only on the linearized 
dynamics (geometrical optics) of the system, one has essentially the same 
stability intervals (in terms of the parameters (s,di,d 2 )) for g.t.p.t.s on the 
surface of constant Gaussian curvature r immersed in the magnetic field /3 
and for g.t.p.t.s on the non-magnetic surface of constant Gaussian curvature 
r e ff = r + (3 2 . As a consequence, one can extend the notions of B-unstable 
and S-unstable g.t.p.t.s to magnetic surfaces of constant curvature as well. 

Equipped with the mentioned above definitions, we are ready now to for- 
mulate the main results of the present paper. Let Q be a billiard table on 
a magnetic surface of constant curvature, and let \(v) > be the Lyapunov 
exponent of the billiard. With any point v = (mi, 9i) 6 V of the phase space 
we associate a formal g.t.p.t. T{y). Let 4>{v) = (m 2 ,9 2 ). We set di = d(v), 
d 2 = d(4>(v)) and s = m\m 2 for the length of particle trajectory between 
mi and m 2 . Then T(v) is determined by the triple (di,d 2 ,s). The formal 
g.t.p.t. T{y) can be realized as an actual g.t.p.t. T(di,d 2 ,s) in an auxil- 
iary billiard table Q v , constructed from the boundary dQ around (see 
||GSG|| ). Let 4>v be the map in Q v corresponding to the g.t.p.t. T(t>), and let 



A(f) = lim^-i-oo - log ||.D0"|| > be the Lyapunov exponent of T(v). Then 
the sufficient condition for hyperbolic dynamics in Q can be formulated as 
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follows (see also Theorem 1 for the alternative formulation in the body of 
the paper). 

Main Theorem. If for fi almost every point v G V, T(v) is B-unstable, 
and for \x almost every point v G V , T(0 n (t>)) is strictly B-unstable for some 
n, then the billiard in Q is hyperbolic (X(v) is positive fi almost everywhere). 



After deriving the conditions, which insure that g.t.p.t.s are B-unstable 
(analogs of (4)), this theorem turns to be a geometric criterion for hyper- 
bolicity of the billiard dynamics. In particular, for planar non-magnetic 
billiards, the Main Theorem yields Wojtkowski's criterion for hyperbolic dy- 
Wo2j |. Let Q be a billiard table satisfying the assumptions of the 



namics 



Main Theorem. Following the approach of Wojtkowski's |[Wo2|| , the metric 
entropy h(Q) = f v X(v) dfi of the billiard can be actually estimated from 
below (Theorem 2 in the body of the paper): 

HQ) > f X(v)dfi. (6) 

The paper is organized as follows. In Section 2 we provide the necessary 
preliminaries and establish the relationship between the geometric optics (i. 
e., the rules of propagation and reflection of infinitesimal light beams) on 
magnetic and non-magnetic surfaces of constant curvature. In Section 3 
we apply these results to obtain explicit analogs of (3-5) for all magnetic 
(non-magnetic) surfaces of constant Gaussian curvature. We obtain a linear 
instability conditions for g.t.p.t.s and show that they distinguish between 
B-unstable and S-unstable trajectories in a natural way. In Section 4, we 
reformulate the Main Theorem in a slightly different way and prove it us- 
ing the invariant cone fields method. We define our cone fields for magnetic 
billiards on all surfaces of constant curvature exactly as in [ |GSG|| . Using 
the geometric optics language, the proof of the preservation of the cone field 
under the assumptions of the Main Theorem is reduced to the corresponding 
non-magnetic problem. The Lyapunov exponent estimation (6) follows from 
the results of ||GSG|| by the same arguments. In Section 5, for each type of 



linearized dynamics we derive the criterion of hyperbolicity for elementary 
billiard tables (the boundary of these billiards consists of arcs of constant 
geodesic curvature). We apply it to construct several classes of magnetic 
billiard tables with hyperbolic dynamics on surfaces of constant curvature. 
Finally, for each type of linearized dynamics we formulate some general prin- 
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ciples for the design of the magnetic billiard tables satisfying the conditions 
of the Main Theorem. Several examples of billiards satisfying these princi- 
ples are also given here. The derivation of the differential condition on the 
boundary of billiards satisfying the Main Theorem is given in the Appendix. 



2 Geometric Optics 

Let M be a surface of constant Gaussian curvature r. We will distinguish 
three cases: M = R 2 -plane (r = 0), M = S 2 -sphere (r > 0) and M = H 2 - 
hyperbolic plane (r < 0). Let us consider the dynamical system arising on M 
from the motion of particle of unit charge, mass and velocity in the presence 
of homogeneous magnetic field perpendicular to the surface. We will denote 
by M(r,/3) the corresponding surface M in the presence of magnetic field 
of strength (3. We can assume, without loss of generality, that (3 > 0. The 
particle trajectories on M(r, /3) are curves of constant geodesic curvature f3 
(circles, paracycles or hypercycles) . We call g s the corresponding flow on 
M(r,(3). 

Let Q be a connected domain in M(r, (3), with a piecewise smooth bound- 
ary dQ. The billiard in Q is the dynamical system arising from the motion 
of a point mass inside Q under the action of magnetic field f3, with specular 
reflections at the boundary. The phase space V of the billiard consists of 
unit tangent vectors, with origin points on dQ, pointing inside Q. The first 
return associated with V is the billiard map, : V — > V. We denote the 
probabilistic invariant measure on V as \x (for its realization in the planar 
case see e.g., ||BK|| ). We will use the standard coordinates (1,9) on V, where 



I is the arclength parameter on dQ and 9, < 9 < n, is the angle between 
the vector and dQ. 

Let D v (p : T V V — > T^^V denote the derivative of <fi. In what follows, we 
are interested in the action of Dcf) on the projectivization B of the tangent 
space TV, see |[Wol|| , ||Wo2|| . The space B = UB V , which is the set of straight 



lines in the tangent planes T V V, v £ V, can be conveniently represented 
using the language of geometric optics. An oriented curve 7 C M, of class 
C 2 , defines a "light beam", i.e., the family of particle trajectories orthogonal 
to 7. The trajectories intersecting 7 infinitesimally close to a point, m £ 7, 
form an "infinitesimal beam" , which is completely determined by the normal 
unit vector n £ T m M to 7, and by the geodesic curvature x °f 7 & t m - Let 
n = v £ V (the point m belongs to the boundary of Q). We denote by B(v, x) 
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the infinitesimal beam determined by the pair (v , x), v G V, x G R- U 00 • 

On the other hand, each pair (t>, x) uniquely defines the line b(v, x) £ B, 
see e.g., ||Wo2|| . Here the vector v defines the corresponding plane T V V, and 



the curvature x defines the direction of the line. Thus, infinitesimal beams 
yield a geometric representation of B = UB V , where x can be used as a 
projective coordinate on the space B v , see [|Wo2|1, ||GSG||. As a result, one 



can study the action of D(f> on B in terms of the action of (f> on the curvature 
of infinitesimal beams B. 

Let p m : T m M — > T m M, m G dQ be the linear reflection about the tangent 
line to dQ. We will use the same letters, 0, p, and g, for the differentials of 
these mappings. Since the billiard map is the composition: 

<i> = p°g, (7) 

it remains to compute the action of g and p on the curvature of infinitesi- 
mal beams. In other words, we need to know how the geodesic curvature of 
an infinitesimal beam changes: a) along free-flight trajectory on M(r,/3); b) 
under reflection. 

a) Bouncless propagation. Let us consider the change of the geodesic 
curvature x( s ) °f an infinitesimal beam along a particle trajectory on M(r, (3). 
As it has been shown in [ Ta2 |, [|Ta3|| (see also [|K| for the planar case) the 



free-flight evolution of x satisfies the Ricatti equation: 

X = r eff + X 2 , (8) 

where r e ff = r + j3 2 is the effective curvature of M(r,/3). For convenience 
we introduce also the parameter £ = |r e jj|^ (then £ -1 is the radius of the 
non- magnetic surface, which has the same linearized dynamics as M(r,/3)). 
Let x — x(0) be the geodesic curvature of an infinitesimal beam at the initial 
point, set x' — x( s ) — 9 s ' X be the geodesic curvature after the time s of 
free-flight evolution. Using eq. (8) one can immediately obtain the relation 
between x an d x' ■ We will distinguish three cases: 

(E) -Elliptic (r eff > 0). M is either S 2 , or R 2 , or H 2 with (3 2 > \r\ 
(strong field). 

x7e = -oot(fro+ f^y s \/ i (9) 

cot(£s) - xl £ 
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(P) -Parabolic (r eff = 0). Mis H 2 and f3 2 = \r\, or M is R 2 and (3 = 0. 



-2 



x=s- 1 + — ; (10) 

s 1 - X 

(H) -Hyperbolic {r e ff < 0). M is H 2 and j3 2 < \r\ (weak field). 

, ,+ i ,j. x sinh~ 2 (£s) . . 

x'n = - coth(e s ) + ^ (ii) 

coth(£s) - xj £ 

Note, for M = H 2 in the case (E) the particle trajectories are circles (tra- 
jectories of finite length), in the case (P) they are paracycles (trajectories 
of infinite length, which touch the boundary of Poincare disc) and in the 
case (H) they are hypercycles (trajectories of infinite length, which have two 
points on the boundary of Poincare disc). 

b) Reflection. Let X-i v - £ T m M be the geodesic curvature and the 
direction of the infinitesimal beam just before the reflection at the point m G 
M. We denote x+ — P • X-i v = v + = P ' v - G ^ to be the geodesic curvature 
and the direction of the infinitesimal beam immediately after reflection. Let 
k be the curvature of dQ at m, and let 9 be the angle between v and the 
positive tangent vector to dQ at m. Then, the extension of the well known 
formula for planar billiards (see e.g., to magnetic billiards on arbitrary 
surfaces (see [ [Ta2| , ||Ta3|| ) gives 



X+ = X-+2K(v), where K(v) = - + . (12) 

smu 

Using classical formulas for surfaces of constant curvature ( |V| , see also [[Tal|| , 
||l'a2j| ), it is possible to give a geometric interpretation of the function K(-). 

Let v E V, and let I G dQ be the origin point of v. Let C{1) C M be the 
osculating circle (resp. paracycle or hypercycle if M = H 2 and \k(1)\ < £) 
of dQ. The free-flight particle trajectory, G(v), corresponding to v intersects 
C(l) at / and at another point I'. In the cases (E) and (P) {r e ff > 0) let d(v) 
be one half of the signed distance between I and I', along G(v). To eliminate 
the ambiguity (when K(v) = 0), we choose the following intervals for d(v): 
d{v) G [—00,00) in the case (P) and £d(v) G [— 7r/2,7r/2) in the case (E). In 
the case (H) (r e // < 0) we will use the following classification of points of 
the phase space V, see fig. 3. We say that v G V is of type A (resp. B) if 
1-^(^)1 > £ (resp. < £). Let V A ,V B be the corresponding subsets of 
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V. Then V = V A U V B is a partition. We denote by d A (v) {d B {v)) one half 
of the signed distance between / and I' along G(v) if v G V A (resp. v G V B ). 
To unite both cases we will use the notation: 



d[v) ~ \d B (v) + q ifveV B . 



Then we have 

K~\v) 



' d(v) in the case (P) 

f-Han^u)) in the case (E) (14) 
_ tanh(£<i(t>)) in the case (H). 



As one can see from eqs. (9-14), the geometric optics (described in terms of 
parameters d, s) depend only on the value r e jj. This fact allows to study 
the linearized dynamics problems on M(r, (3) using the corresponding results 
for the non-magnetic surfaces of the constant Gaussian curvature r e ff = 
r + [3 2 . The corresponding transition M(r,/3) M(r e jj,0) is schematically 
illustrated by fig. 4. 



3 Stability of Generalized Two- Periodic Tra- 
jectories 

Let Q be a billiard table on M(r,(3). For each v G V let t(v) be the corre- 
sponding past semitrajectory in Q. Consider the curvature evolution of an 
infinitesimal beam along t(v). Starting with B(<f)~ k ■ v, x) for arbitrary x, we 
obtain after k steps forward the infinitesimal beam B(v,x^), X^ = 4> k ' X- 
Eqs. (9-11) and (12) describe the action of the billiard map on the curva- 
ture of infinitesimal beams. Assuming k to be infinity, we obtain a formal 
continued fraction corresponding to the semitrajectory t(v): 

c(v) = X {oo) = a + — h . (15) 

0-1 

a_i H 

a_ 2 • • • 

The coefficients of the continued fraction are determined by di = d{$ % ■ v), 
and by the lengths of consecutive billiard segments as follows: 
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(P) a i = -2sj 1 + 2dj\ k = -sr 2 ; 

(E) cii = -2 cot(£sj) + 2 cot(^), bi = - sin~ 2 (£,s;); 

(H) di = -2 coth(^Si) + 2 coth(^dj), h = - sinh~ 2 (^j). 

The continued fractions (15) determines the stability type of the trajectory: 
t{v ) is unstable if c{v ) is convergent (see e.g., [EI). Since for a given sequence 
of di and Sj, c{v ) is completely determined by r e ff, one can reduce the prob- 
lem of stability trajectories on M(r,/3) to the corresponding "non-magnetic" 
problem on M(r e ff, 0). 

As it has been mentioned in the introduction, we are interested in the 
stability properties of generalized two periodic trajectories. A trajectory is a 
generalized two periodic trajectory (g.t.p.t.) if its parameters di are periodic: 
<i 2 i+i = di, d 2i = d 2 and Si = s are the same along the trajectory (see fig. 
2). Obviously, a g.t.p.t. yields a periodic continued fraction. We denote 
by T(di,d2,s) the g.t.p.t. with parameters (di,d 2 , s) and by c{di,d 2 ,s) the 
associated continued fraction. 

The stability of T(di, d 2 , s), or equivalently, the convergence of the two 
periodic continued fraction c(<ii, <i 2 , s) has been studied in ||GSG|| for non- 



magnetic surfaces of constant curvature. On the basis of the equivalence 
between the magnetic and non-magnetic problems we can immediately gen- 



eralize the results of ||GSG|| to the case {3 ^ 0. 

Proposition 1. The continued fraction c(di,d2,s) converges if and only if 
the following inequalities are satisfied. 

(P) (s - di) (s - d 2 ) (s-di- d 2 )s > 0; 

(E) sin(£(s - di)) sin(£(s - d 2 )) sin(£(s - d x - d 2 )) sin(£s) > 0; 

(H) sinh(£(s - dj) sinh(£(s - d 2 )) sinh(£(s - d 1 - d 2 )) sinh(£s) > 0. 

Below we reformulate Proposition 1 explicitly as conditions for the insta- 
bility of the corresponding g.t.p.t. 
(P) T(di, d 2 , s) is unstable if and only if 

![d±, d 2 ] U [d\ + d 2 , oo) if d%, d 2 > 
[0,oo) ifd 1 ,d 2 <0 (16) 

[0, d x + d 2 ) U [d u oo) if d 1 > 0, d 2 < 0. 

(E) In this case > £s > 2tt, and we set n = n ■ 
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Then T(di, d 2 , s) is unstable if and only if 

{[di + d 2 , ft] U [di, d 2 ] if di, d 2 > 

[0, d\ + d 2 + 7f ] U [ft — d\ , 7f — d 2 ] if di, d 2 < 
[d 2 , * + di] U [0, di + d 2 ] if di < 0,d 2 > 0, \d 2 \ > \d x \ 

[d 2 , ft + di] U [ft + d 2 + c?i, 7f] if di < 0, d 2 > 0, | c?2 1 — \di\- 

(17) 

(H) It matters whether Vi G V A or t>j G for i = 1,2. We say that 
T(di,d 2 ,s) is of type (A — A) if Vi G V" 4 and t> 2 G VA The other types: 
(A — B), (B — A), and (B — B) are defined analogously. We formulate the 
explicit criteria of instability for T(d±, d 2 , s) type-by-type. 
Type (A — A): 

U [df + df, oo) if df, df > 

if<4<0 (18) 
[0, df + df] U [df, oo) if df >0,df< 0, 




s£ i ri T« 2 )^ if df + df > , x 

l rn if df + df < 0, 1 ; 




Type (5 - B) 



Types (A - B) or (B - 

_f[#,oo) if ^ > 

G \[0,oo) if#<0, ^ 

It is worth mentioning that in Proposition 1 (resp. eqs. (16-20)) the hyper- 
bolicity of T(d 1 ,d 2 ,s) corresponds to strict inequalities (resp. inclusions in 
the interior). The equality case (resp. boundary case) corresponds to the 
parabolicity of T(d\, d 2 ,s). There are also two special cases when T(di, d 2 , s) 
is parabolic independently of the value of s: (P), d\ = d 2 = — oo and (H), 
\df\ = \df\ = oo. 

We call the right hand side of eqs. (16-20) the instability set of T(di, d 2 ,s). 
In general, it is a union of two intervals, where one of them degenerates when 
|di| = |d 2 |, while the other is always nontrivial. Following the terminology 
of our previous work IpSQI , we will say that the interval which persists 
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is a "big interval", while the other one is a "small interval". We will say 
that T(di,d 2 ,s) is (strictly) B-unstable if s belongs to the (interior of the) 
big interval of instability. The proposition below makes this terminology 
explicit. 

Proposition 2. The g.t.p.t. T(di,d 2 ,s) is B-unstable if (and only if) the 
triple (di,d 2 ,s) satisfies the following conditions: 

(P) 

[di + d 2 , oo) if di, d 2 > 

s e { [0,oo) if di,d 2 < (21) 

[di,oo) ifdi>0,d 2 <0 



(E) 



[d% + d 2 , it] if di, d 2 > 
s modvf e { [0, di + d 2 + 7f] ifdi,d 2 <0 (22) 
[d 2 , 7T + dj if di < 0, d 2 > 



(H) The case (A - A) 



[df + df, oo) if df, d£>0 
s G | [0, oo) if df , d^ < (23) 

[df , oo) if df > 0, do 4 < 0, 



or |d^| = | d 2 \ = oo and arbitrary s. 
(H) The case (B - B) 



f[df + d?,oo) ifdf + df>0 

fc \[0,oo) ifdf + df<0 1 j 

(H) The cases (A - B) or (B - A) 

/K\oo) if ^ > ( , 

SG \[0,oo) ifd^<0. [Zb) 

Obviously, the conditions (21-25) for B-unstable g.t.p.t.s are the same as 
those which appeared in for the corresponding non-magnetic cases. 



4 The Main Theorem 

Let Q be a billiard table and v G V be an arbitrary point in the phase space 
of the billiard map. Set V\ = v , v 2 = (j)(v), di = d(v j), i — 1, 2, and let s = s(t>) 
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be the length of the particle trajectory between the origin points of V\ and v 2 
respectively. We will associate with v a formal g.t.p.t. T(v) = T(di, d 2 , s), 
which parameters are defined by the triplet (d\, d,2, s). We denote by X(v) 
the Lyapunov exponent of the billiard Q and by A(t> ) the Lyapunov exponent 
of T{v ) (see Sect. 1), which are defined for /^-almost all v G V. 

Using Proposition 2 we introduce the following special class of points of 
the phase space of the billiard map. 



Definition 1. A point v G V of the billiard phase space is 

a) B- hyperbolic (or strictly B-unstable) if the corresponding g.t.p.t. T(v) 
is strictly B-unstable; 

b) B-parabolic if the corresponding g.t.p.t. T(v) is B-unstable and parabolic 
(i.e., s belongs to the boundary of the appropriate interval (21-25) ); 

c) B-unstable if the corresponding g.t.p.t. T(v) is B-unstable (i.e., B- 
parabolic or B-hyperbolic) ; 

d) eventually strictly B-unstable if there is some integer n such that 
T((p l (v)) is B-unstable for < i < n and T(<f) n (v)) is strictly B-unstable. 



Below we formulate the main theorem of the present work. 



Theorem 1. Let Q be a billiard table on M(r,/3). If fi- almost every point of 
the billiard phase space is eventually strictly B-unstable, then the Lyapunov 
exponent A is positive ^-almost everywhere. 



Proof. The proof of the theorem is based on the cone field method which has 
been initially applied to the planar billiards in || Wo 1|| 



]Wo2 | 



A cone in T V V corresponds to an interval in the projectivization B v . 
Therefore, a cone field, W, is determined by a function, W(-), on V, where 
each W(v ) is an interval in the projective coordinate x- We define the func- 
tion W{v ) as in ||GSG|| . For completeness, we repeat this definition below. 



(P) and (E) 



(H) 



W(v) 



W(v) 



([K(v),+oo] 
\[-oo,K(v)] 

([K{v),+oo] 
\[-oo,K(v)} 



if K(v) > 
if K{v) < 



if K{v) 
if K(v) 



As it follows from Lemma 2 in ||GSG|| , this cone field is eventually strictly 
preserved by the billiard map if the conditions of Theorem 1 are satisfied. 



14 



By this fact the proof of the theorem follows immediately from Wojtkowski's 
theorem (Theorem 1 in |[Wo2|| ) . □ 



Applying the method developed in [ |Wo2j| , one can actually estimate from 
below the Lyapunov exponent using the cone field defined above. The result 
is given by the next theorem. 

Theorem 2. Let Q be a billiard table satisfying the assumptions of Theorem 
1, then 

h{4>) = / Mv) dfx> X(v) d/i. 
Jv Jv 



Proof. The proof follows immediately by the repetition of calculations given 
in the proof of the analogous theorem for the non-magnetic case (see Theorem 
2 in flGSGfl ). □ 



5 Applications and Examples 

Theorem 1 together with Proposition 2 lead to a simple geometric criterion 
for billiard tables with hyperbolic dynamics. In this section we apply this 
criterion to construct various classes of hyperbolic billiards on M(r, f3). 

5.1 Elementary billiard tables 

There is a class of billiard tables, where the application of Theorem 1 gives 
an especially simple criterion for hyperbolicity. This class consists of billiard 
tables Q, whose boundary is a finite union of arcs, Tj, of constant geodesic 
curvature, K(Ti) = Ki. We call these tables elementary. We will use the 
notation Tf (resp. T~) if K(Ti) > (resp. n(Ti) < 0). Let Ci be the 
curve of constant geodesic curvature such that Tj C d and Di C M be 
the corresponding disk {Ci = dDi). Since the representation dQ = U^ 1 Tj 
is unique, we call Tj the components of dQ. In the following, we consider 
elementary billiard tables for which > (3. One may easily see that the 
fulfillment of this inequality is necessary for billiards satisfying the conditions 
of Theorem 1 (see discussion in the Section 5.2 for billiards with boundaries 
of general type). 
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(E) Elliptic case (r e ff > 0). Let D C M be a disc such that <9.D is the 
circle whose geodesic curvature k satisfies k > (3. We define the component 
— D C M as set of the points m! G M satisfying the condition m/m = 7f 
for some point m G -D, where m'm is the length of the particle trajectory 
between the points m, m'. We will refer to — D as the dual component of 
D = +D. Straightforward analysis shows that — D is the ring whose width 
equals to the diameter of D and its radius is defined by £ (for M = R 2 its 
radius is see figs. 5a,b,c. When M = S 2 and (3 = 0, —D is the disk 

obtained from D by reflection about the center of S 2 , as it has been defined 
in pgq . 

Let us also introduce the terminology: If R C S C M are regions with 
piecewise C 1 boundaries, we call an inclusion R C S proper if dRHint S ^ 0. 
The application of Theorem 1 to the elementary billiard tables in the case 
r e ff > leads to the following criterion for hyperbolicity. 

Corollary E. Let Q C M be an elementary billiard table whose boundary 
consists of N > 1 components of type plus or minus. Suppose Q satisfies the 
following conditions: 

Condition El. For every component Tf of dQ we have A C Q. Besides, 
either —Di C Q, or —Di C M \ Q, where the inclusions are proper; 

Condition E2. For every component Tj we have Dj C M \ Q, and the 
inclusions —Dj C M\Q, or —Dj C Q are proper. 
Then the billiard in Q is hyperbolic. 

Outline of proof: The assumptions of Corollary E imply those of Theorem 1. 
Remark. Suppose Q' = M \ Q is connected. If Q satisfies Conditions El and 
E2, then Q' also does, and hence the billiard in Q' is hyperbolic. 

Examples. "Lorenz gas" billiards. Such billiards are obtained by removing 
from M a number of disjoint discs Di, so that Q = M \ UDj. If all the 
intersections Di D ±Dj % ^ j, are empty, then the billiard in Q is hyperbolic 
by Corollary E. The simplest example of such hyperbolic billiard is obtained 
by removing two disks from the magnetic plane, see fig. 6a. 

The intersections Di n ±Dj % ^ j, are always empty, if all the discs are 
contained inside of a free-flight particle trajectory (i.e., if all the discs lie 
inside a circle of geodesic curvature (3). Such billiards are the "magnetic" 
analogs of the non-magnetic hyperbolic billiard tables on the sphere, obtained 
by removing a finite number of disjoint disks from one hemisphere ||GSG|| . 
The examples of hyperbolic billiards of this type on S 2 , R 2 and H 2 are shown 
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in fig. 6b,c,d. 

One can consider also unbounded billiard tables Q obtained by removing 
an infinite number of disjoint disks from R 2 , H 2 . The simplest example of 
this type is obtained by removing a chain of equal disks from M = R 2 , as 
shown in fig. 7a (this billiard can be also seen as cylinder with one hole). 
Because of the translation symmetry, one needs to check the non-intersection 
condition only for one disk. The non-intersection condition is also necessary 
for hyperbolicity of such billiards. If it is not satisfied, then Q has at least 
two stable g.t.p.t.s (see fig. 7a). 

Another type of unbounded hyperbolic billiard tables can be obtained 
by removing a lattice of the disks from M = R 2 , H 2 . The example of such 
billiard shown in fig. 7b, is equivalent to the torus with one hole. Here, again, 
because of the translation symmetry, one has to check the non-intersection 
condition only for a single disk. 

"Flowers" like billiards. Consider a simply connected billiard table Q, 
whose boundary consists of several circular arcs of positive and negative 
curvature satisfying the condition |/^| > (3. Such billiards were originally 



introduced by Bunimovich [Bul|| Pu2|| as examples of planar (non-magnetic) 



hyperbolic billiards with convex boundary. It has been demonstrated that 
for r = 0, (3 = such billiards are hyperbolic if the conditions Di C Q are 
satisfied for each convex component of the boundary. For r e ff > we have by 
Corollary E the additional requirement: dQ D —Di = for each component 
of the boundary (compare with the analogous conditions in [GfSG] for the 



case (3 = 0, M = S 2 ). The examples of hyperbolic billiards Q on S 2 , R 2 , 
H 2 of "flower" type satisfying the conditions of Corollary E are shown in 
fig. 8a,b,c. It follows from the remark above that billiards in the domain 
Q' = M\Q are also hyperbolic. 

(H+P) Hyperbolic and parabolic cases {r e ff < 0). The criterion for 
hyperbolicity in this case is given by the following corollary. 

Corollary H. Let Q C M be an elementary billiard table, and let dQ consist 
of N > 1 components. If Q satisfies conditions: 

Condition HI. For every convex component Tf of dQ, we have Di C Q; 

Condition H2. For every concave component of dQ, we have ft(r~) < —f3 
and for every convex component n{Tf) > £. 
Then the billiard in Q is hyperbolic. 

Outline of proof: The assumptions of Corollary H imply those of the Theorem 
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1. 

Remark. When (3 — > and r — > 0, the condition H2 is automatically fulfilled 
and Corollary H turns to be the classical criterion of Bunimovich ||Bu2|| for 
hyperbolicity of planar, non-magnetic billiard tables. 



Examples. Analogs of Sinai billiards. The boundary of these billiards con- 
sists of concave arcs T~ of constant curvature (see fig. 9). If the condition 
ft(IT) < —j3 is satisfied for each component of the boundary, then the billiard 
is hyperbolic by Corollary H. 

Analogs of Bunimovich billiards. The example of hyperbolic billiard table 
with convex components satisfying the conditions HI, H2 is shown in fig. 10. 
Remark. The assumptions in Corollaries E and H that N > 1 and that the 
inclusions be proper are needed only to exclude certain degenerate situations, 
where each v G V is B-parabolic. This is the case, for instance, if Q is a disc, 
or the annulus between concentric circles. 



5.2 Hyperbolic billiard tables with boundary of gen- 
eral type 

Let us consider billiard tables on M(r,/3) with piecewise smooth boundary, 
dQ = Uj7i of general type. The components 7j are C 2 smooth curves parame- 
terized by the arclength I, whose curvature Ki(l) has the same sign along each 
7i. We will refer to 73 as convex component if K{(1) > 0, or as concave com- 
ponent if Ki(l) < 0. Let us denote «(7i) = max{Kj(7),Z G 7^} for the convex 
components, and K(ji) = min{Kj(7),Z G 7^} for the concave components. 

Following the terminology in |[Wo2|| , we introduce the class of convex 
scattering curves on M(r,/3). 

Definition 2. A smooth convex curve 7 C M is (strictly) convex scattering 
if for any v G V, such that the origin points of v and <p(v) belong to 7, the 
corresponding g.t.p.t. T(v) is (strictly) B-unstable. 

A curve 7 is convex scattering if one of the relevant conditions (21-25) 
is satisfied for each pair of points on 7. Regarding the planar non-magnetic 
case, this leads to the definition of Wojtkowski [|Wo2j| for convex scattering 
curve. Let us introduce the parameter R(l) = (k(1) — Considering 
the infinitesimally close points on 7 we show in Appendix that the condition 
R"(l) < is necessary for 7 to be convex scattering. It should be noted, that 
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this condition is also sufficient in the planar, non- magnetic case (see ||Wo2|| ), 
but not for generic parameters r, (3 (see |pSG|| for f3 = case). 

In what follows, we formulate the principles for design of hyperbolic bil- 
liards satisfying the conditions of Theorem 1. Let Q be a billiard table 
satisfying the conditions of Theorem 1. Then each convex component of dQ 
has to be convex scattering and consequently, the condition R" < holds 
along each convex component of the boundary. There is an additional re- 
striction on the curves 7« which compose the boundary of Q. It follows from 
Proposition 2 that for billiards satisfying the conditions of Theorem 1 the 
sign of K(v) (d{y)) depends only on the origin point of v (there is no depen- 
dence on 9) for any v G V, i.e., K(v) (d(v)) has the same sign along 73 as 
re(7i). This happens if for each component ji, |re(7i)| > (3 (the magnetic field 
is sufficiently weak). Thus, in what follows we particularly exclude from our 
consideration the magnetic billiards with flat boundaries. Such billiards do 
not satisfy the conditions of Theorem 1. 

Design of hyperbolic billiard tables in the (E) case. 

By Definition 2 a curve 7 is convex scattering if it is convex and the condition 

di + d 2 < s < 7T, (26) 

holds for any pair of points on 7. For simplicity of exposition, we will restrict 
our attention for M = R 2 , H 2 to the bounded billiard tables and for M = 
S 2 to the billiard tables which can be placed in a hemisphere. Theorem 1 
yields the following principles for the design of piecewise billiard tables with 
hyperbolic dynamics in (E) case: 
PI: |k(7i)| > (3 for all components. 

P2: All convex components of dQ are convex scattering. 
P3: Any convex component of dQ has to be "sufficiently far", but not "too 
far", from any other component. Any concave component has to be not "too 
far", from any other concave component. 

The precise meaning of P3 is that the parameters of any two consecu- 
tive bouncing points, which belong to different components of the boundary, 
satisfy the condition (22). In particular it implies the set of restrictions 
on the angles between consecutive components of the boundary. It can be 
formulated as an additional principle. 

P4: Let 7i,7i+i C dQ be two adjacent components, meeting at a vertex. If 
both 7j and are convex, then the interior angle at the vertex is greater 
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than 7T. If 7j and 7j + i have different sign of curvature, then the angle in 
question is greater or equal to ir. 

Another restriction which arises from P3 is that the length (equivalently 
the time) of free-flight between any two consequent bouncing points on the 
boundary of the billiard has to be not greater than n. In other words, the 
billiard table has to be "smaller" than circle drawn by a free-flight particle 
on M(r,ft). 

Examples. The examples of the hyperbolic billiards on R 2 satisfying the 
above principles are shown in fig. lla,b. A bounded Sinai-like billiard, whose 
boundary consists of (strictly) concave components (fig. 11a) always satisfies 
the principles P1-P4 for sufficiently weak magnetic field. 

The example of a convex billiard is shown in fig. lib. It is a cardioid, 
whose boundary is strictly convex scattering curve for ft = (see [ Wo2[j ). For 



ft = this billiard is hyperbolic, as it follows from Theorem 1. Since strictly 
convex scattering curve remains to be such under small perturbations of /3, 
the billiard in fig. lib is hyperbolic for sufficiently weak magnetic field. 

Design of hyperbolic billiard tables in the (P+H) case. 

Definition 2 leads to the following geometric conditions on the convex scat- 
tering curve in the (P+H) case. A convex curve 7 is convex scattering if 
« (7) > £ and for each pair of points on 7 

d x + d 2 < s. (27) 

Theorem 1 yields the following principles for the design of billiard tables 
with hyperbolic dynamics in the (P+H) case: 

PI: «(7i) > £ for any convex component of dQ and K(ji) < —f3 for any 
concave component of dQ. 

P2: All convex components of dQ are convex scattering. 
P3: Any convex component of dQ is "sufficiently far" from any other com- 
ponent. 

More precisely, condition P3 means that any two consecutive bouncing 
points of the billiard ball, which belong to different components, satisfy eqs. 
(23-25). In particular, this yields, the same inequalities (P4) as in (E) case, 
for the interior angles between consecutive components of dQ. 

Examples. In (P+H) case, any concave billiard is hyperbolic if the condition 
< — ft is fulfilled for each component of the boundary. As in the case 
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(E), the examples of the convex hyperbolic billiards can be obtained from 
their non-magnetic counterparts satisfying the conditions of Theorem 1. 

Finally, it should be noted, that formulated above principles for design 
of hyperbolic billiards on M(r,/3) are robust under small perturbations of 
j3, r and the billiard wall. Generally, one can construct hyperbolic billiards 
on magnetic surfaces of constant curvature on the basis of the corresponding 
non-magnetic planar billiards satisfying Wojtkowski's criterion. 



6 Conclusions 

In the present paper we have formulated the criterion for hyperbolic dynam- 
ics in billiards on surfaces of constant Gaussian curvature r in the presence of 
a homogeneous magnetic field /3 perpendicular to the surface. The criterion 
is valid for all values of r, f3 and its geometric realization depends only on the 
type of linearized dynamics (elliptic, parabolic or hyperbolic). In this way we 
extend our recent results in ||GSG|| to the case of magnetic surfaces of constant 
curvature. The basic property, which allows unification of the hyperbolicity 
criteria for the magnetic and non-magnetic billiards on surfaces of constant 
curvature, is the equivalence between the geometric optics in both cases. In 
fact, in terms of special parameters di, s, the geometric optics depend only on 
the effective curvature r e ff = r + j3 2 of the surface. It is important to stress, 
that the dynamics in magnetic and non-magnetic billiards are very different 
(e.g., the magnetic field breaks time reversal symmetry). It is the only lin- 
earized dynamics, which are the same for the considered systems. Applying 
the hyperbolicity criterion, we were able to construct the different classes of 
hyperbolic billiards for each type of the linearized dynamics (equivalently for 
each of the signs of r e ff). 

There are two types of necessary conditions which arise for hyperbolic 
billiards satisfying our criterion. The first one is a requirement for the convex 
components of the boundary to be convex scattering. As a consequence, the 
inequality R"(l) < has to be satisfied along each convex component. This 
inequality is generalization of well-known Wojtkowski's condition [ Wo2 | for 



convex component of planar (non-magnetic) hyperbolic billiard. It has been 
demonstrated for planar non- magnetic billiards in jBu3|, |[Bu4 | , [ Do] that 



Wojtkowski's criterion can be considerably strengthened. This suggests, in 
particular, that condition R"(l) < can be relaxed for general parameters r, 
(3 by employing invariant cone fields, different from the one used in the present 
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paper (see discussion in ||GSG|| ). The second type of conditions is specific for 
magnetic billiards. This is a requirement of "weakness" for the magnetic field 
compared to the curvature of the billiard boundary. For generic systems, 
such condition is expected, in order to prevent stable skipping orbits close 
to the boundary. It has been shown in [BR] , (see also |BK| ) that billiard 



with sufficiently smooth boundary possesses invariant tori corresponding to 
skipping trajectories. It seems that in the strong field regime a part of 
stable periodic orbits has to survive even if the smoothness of the boundary 
is broken. It remains, however an open question, whether the condition 
\Ki\ > (3 can be relaxed for generic billiard. 

The positive Lyapunov exponent for a billiard implies strong mixing prop- 
erties: countable number of ergodic components, positive entropy, Bernoulli 
property etc. It should be pointed out, however, that ergodicity does not 
automatically follow from the positivety of Lyapunov exponent. Neverthe- 
less, one can expect that billiards satisfying the conditions of Theorem 1 
will be typically ergodic. It seems that the methods developed for the proof 
of ergodicity of planar hyperbolic billiards can be extended to the class of 
billiards considered in the present paper. 
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7 Appendix 

We will investigate the conditions under which a convex arc on the surface of 
constant curvature M in the presence of magnetic field (3 is convex scattering. 

For simplicity of exposition, we consider the case, when M is magnetic 
plane. Let 7(Z) C M be any smooth curve parameterized by arclength /, 
and let k(Z) be the geodesic curvature of 7. Let now 7(^0) an d 7(^1) be two 
points on 7, such that the arc of 7 between 7(^0) and 7(^1) lies entirely on 
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one side of straight line passing through 7(^0) an d 7(^1)- We choose cartezian 
coordinate system (x,y) in such a way that y(lo) = y(h) — 0, x(Iq) = —x(li) 
and the arc of 7 between 7(^0) and 7(^1) lies above a;- axis, see fig. 12. Let 
a(l) be the angle, which ^ makes with x-axis , then 

dx dy . da . . 

— = cos a; — = sin a; — = —k. (28) 
dl dl dl 

We introduce also an auxiliary variable 5, such that fix = sin 5. 

For (3 > there are two different particle trajectories connecting the 

points 7(/ ) and 7(^1) (resp. two different g.t.p.t.s corresponding to these 

points), see fig. 12. Below, we consider the trajectory which lies in the 

lower halfplane. Then, the results for trajectory in the upper halfplane are 

obtained by the change of the sign of f3 to the opposite. Let 9 = a + 5. Then, 

at the points Z ,i> ^0,1) are the angles between 7 and the particle trajectory 

connecting 7(7 ) and 7(^1). Set A = s — di — d 2 - By eq. 14 we get 



A = /T 1 / \d{ arctan ( J^LX) + d5 
J \ \ k — 3 cos 9 



\ \k — (3cos9 / 

r /-K'sm9 + K(K + s -^)(^-cos9)\ 

_ J] I V _ smS J \ cosS ) / hq\ 

~J dl [ k 2 + 3 2 -23kcos9 ) [Zy) 

We separate the last integral into the sum of two parts. The first one is 




-re' sin \ r . / R' 



i 2 + p 2 - 23k cos 9 J J ~* \ 1 + AR 2 k[3 sin 2 9/2 J 
where R~ l (l,f3) = k(1) — 13. Since y(lo) = y(h) = 0, we obtain 

y R'(AR 2 Kp sin 2 §)' \ R"L 3 k , 



+ 4 J R 2 «/3sin 2 § (l + 4i?%3sin 2 f)V 12 + 0(L )? 

(30) 

where L — li — l is the length of the curve between the points 7(^0), 7(^1)- 
Analogously, for second part we have 



// = dl 



'k(k + (3 s -^) fseas - C os^ 

\ " sin J \ cos j 

k 2 + 3 2 - 23kcos9 
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, i cos<5 V sin a ^ i rM T i\ 

dy| ^ + ^-2^oosg l=Q(L } " 
Adding both parts we obtain finally 

A = / + // = -^ + 0(L 4 ). (31) 

Thus, if the curve 7 is convex scattering, then the condition R"(l, /3) < holds 
everywhere on 7. Considering trajectories of the second type (i.e., trajectories 
which lay in the upper halfplane), we obtain the condition R"(l, —(3) < for 
convex scattering curves. However, it easy to see, that R"(l,/3) < actually 
implies R"(l, -/?) < 0. 

Repeating the same analysis for general M(r,fl) we have found (see also 
|GSG| for (3 = case) that eq. 31 holds for all surfaces of constant curvature. 



As a consequence, R" < is a necessary condition for convex scattering on 
M(r,/3). On the contrary, if the strict inequality R" < holds along 7, then 
by eq. 31, any sufficiently small piece of 7 is convex scattering. 
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fig. 11a 



20 




-3* 



fig. lib 



21 



a y 




